Ah initio no-core full configuration calculations of light nuclei 
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We perform no-core full configuration calculations for a set of light nuclei including ^^O with a 
realistic NN interaction, JISP16. We obtain ground state energies and their uncertainties through 
exponential extrapolations that we demonstrate are reliable in ^H, '^H and *He test cases where 
fully converged results are obtained directly. We find that ^He, ®Li and *He are underbound by 
about 600 keV, 560 keV and 1.7 MeV, respectively. ^^C is overbound by about 1.7 MeV and ^'^O 
is overbound by about 16 MeV. The first excited 0"'' states in ^^C and ^^O are also evaluated but 
their uncertainties are significantly larger than the uncertainties for the ground states. 

PACS numbers: 21.60.Dc, 21.60.Cs, 21.45.-v, 21.30.-x, 27.20.-|-n, 27.10.+h, 21.10.-k, 21.10.Dr 



I. INTRODUCTION AND MOTIVATION 

The rapid development of ab initio methods for solv- 
ing finite nuclei has opened a range of nuclear phenom- 
ena that can be evaluated to high precision using realistic 
nucleon-nucleon (NN) and three-nucleon (NNN) interac- 
tions. Such advances define a path for testing fundamen- 
tal properties of the strong interaction such as their ori- 
gins from QCD via chiral effective field theory [1-4]. In 
addition, they prepare a foundation for nuclear reaction 
theory with unprecedented predictive power. 

Here we investigate the direct solution of the nuclear 
many-body problem by diagonalization in a sufficiently 
large basis space that converged binding energies are ac- 
cessed — either directly or by simple extrapolation. Our 
choice is a traditional harmonic oscillator (HO) basis so 
there are two basis space parameters, the HO energy hfl 
and the many-body basis space cutoff Nmax ■ Nmax is de- 
fined as the maximum number of total oscillator quanta 
allowed in the many-body basis space above the mini- 
mum for that nucleus. We obtain convergence in this 
2-dimensional parameter space {Ml,N„iax), where con- 
vergence is defined as independence of both parameters 
within evaluated uncertainties. 

Since we treat all nucleons equivalently and we achieve 
convergence within evaluated uncertainties, we refer 
to our approach as the "No-Core Full Configuration" 
(NCFC) method. The NCFC is both related to and dis- 
tinct from the No-Core Shell Model (NCSM) [5], that fea- 
tures a finite matrix truncation, and an effective Hamil- 
tonian renormalized to that finite space. The regulator, 
Nmax, appears in our NCFC, where it is taken to infinity, 
and in the NCSM, where it also appears in the definition 
of the effective Hamiltonian. In both NCFC and NCSM, 
this choice of many-body basis regulator, N^ax, is needed 
to preserve Galilean invariance — to factorize all solu- 
tions into a product of intrinsic and center-of-mass mo- 
tion components. With Nmax as the regulator, both the 
NCFC and the NCSM are distinguished from the "Full 
Configuration Interaction" (FCI) method in atomic and 
molecular physics that employs a cutoff in single-particle- 
space. 



The NCFC results should agree with the NCSM and 
no-core FCI results when the latter results are obtained 
in sufficiently large basis spaces. In the case of NCSM, 
larger cluster sizes for the effective Hamiltonian may be 
employed to accelerate convergence. 

Given the rapid advances in algorithms and comput- 
ers, as well as the development of realistic non-local NN 
interactions that facilitate convergence, we are able to 
achieve converged results, either directly or through ex- 
trapolation, without recourse to rcnormalization of the 
interaction. That is, with our adopted interaction, we 
do not need to soften the NN interaction by treating it 
with an effective interaction formalism. Rcnormalization 
formalisms necessarily generate many-body interactions 
that significantly complicate the calculations and are of- 
ten truncated for that reason. Renormalization without 
retaining the effective many-body potentials abandons 
the variational upper bound characteristic that we pre- 
fer to retain. Furthermore, convergence with increasing 
model space is generally neither uniform nor monotonic 
when applying renormalization without retaining the in- 
duced many-body potentials. This leads to challenges for 
extrapolation to infinite model spaces. 

Our NCFC approach requires methods to reliably ex- 
trapolate results obtained in a finite basis space to the 
infinite or complete basis space limit. This need for ex- 
trapolation tools mirrors similar situations in other fields 
of science where a sequence of results with increasing 
resolution must be extrapolated to the limit of infinite 
resolution. The resulting high precision results of the 
NCFC do not agree exactly with experiment. Indeed, no 
known realistic NN interactions provide exact descrip- 
tions of a similar range of nuclear data that we examine 
and it is probable that NNN and higher body interactions 
are needed. 

In order to further motivate our efforts to develop ro- 
bust extrapolation tools, we show in Fig. 1 the Hamilto- 
nian matrix dimensions for a set of representative light 
nuclei. We employ the "m-scheme" where each HO 
single-particle state has its orbital and spin angular mo- 
menta coupled to good total angular momentum, j, and 
magnetic projection, m. The many-body basis states are 
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FIG. 1: (Color online) Representative Hamiltonian matrix di- 
mensions for total magnetic projection M = states in the 
single-particle m-scheme as a function of the maximum to- 
tal oscillator quanta of excitations specified by Nmax- The 
natural parity matrix dimensions are represented by the spe- 
cific points while unnatural parity matrix dimensions would 
he close to the interpolating lines at odd values of Nmax- 



Slater determinants in this HO basis and are limited by 
the imposed symmetries — parity and total angular mo- 
mentum projection (Af), as well as by Nmax- In the 
natural parity cases for even nuclei shown, M ~ en- 
ables the simultaneous calculation of the entire spectrum 
for that parity and Nmax- 

The nearly exponential growth in matrix dimension 
with increasing Nmax is clearly evident in Fig. 1. In 
order to achieve NCFC results for the heavier of these 
nuclei by extrapolation, using realistic interactions, we 
would need to diagonalize matrices that are beyond the 
reach of present technologies. However, in cases up to and 
including ^^O, we may expect to obtain systematic results 
for the first few increments of Nmax- In order to utilize 
the sequence of results obtained with Nmax values that 
are currently accessible, we need to investigate suitable 
extrapolation tools. 

To better understand the scale of computational effort 
needed for no-core microscopic nuclear structure studies, 
we consider the memory storage demands as a function 
of matrix dimension. For several representative nuclei, 
we enumerate the number of non-vanishing matrix ele- 
ments of the resulting many-body Hamiltonian matrix 
(its lower triangle only for efficiency) and display the re- 
sulting counts as a function of the matrix dimension in 
Fig. 2. We present results for the case of a 2-body input 
Hamiltonian (NN interaction only) and for the case of a 3- 
body Hamiltonian (NN -I- NNN interactions). In spite of 
the very large memory requirements, the various curves 
display an encouraging trend. Specifically, the number 
of non-vanishing many-body matrix elements follows a 
£)3/2 growth rate, where D is the dimension of the ma- 
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FIG. 2: (Color online) Number of non-vanishing many-body 
Hamiltonian matrix elements for representative light nuclei 
as a function of the basis space dimension. The points rep- 
resent sample cases that have been solved and correspond to 
those indicated in the legend. The curves approximate a 

^3/2 

power law where D is the basis space dimension. The vertical 
arrow measures a factor of 30 between the 2-body Hamilto- 
nian (rank=2) and 3-body Hamiltonian (rank=3) cases for 
^^C at the same dimension corresponding to Nmax = 6. Note 
the logarithmic scales. 



trix. That is, the matrices exhibit a very sparse character 
and this is the property that allows us to diagonalize the 
large matrices that we can presently solve. 



II. SELECTION OF HAMILTONIAN 
INGREDIENTS AND EXAMPLE NUCLEI 



In order to carry out the NCFC calculations, we 
require a realistic NN interaction that is sufRciently 
weak at high momentum transfers that we can obtain 
a reasonable convergence trend. The conventional Lee- 
Suzuki-Okamoto renormalization procedure of the ab 
initio NCSM [5] develops soft, Nmax-dependent, effective 
interactions that provide answers close to experimental 
observations. However, the convergence trend of the re- 
sults with increasing Nmax is often not uniform and leads 
to challenges for extrapolation to infinite model spaces. 
Nevertheless, there is also encouraging progress in ex- 
trapolating NCSM ground state energies of light nuclei 
using different strategies [6, 7]. Of course, as the basis 
space increases, one expects the NCSM and NCFC meth- 
ods to arrive at the same exact result. Thus, the choice 
of method, NCSM or NCFC, will ultimately depend on 
the underlying Hamiltonian selected for the application. 
In the NCFC approach discussed here, we seek to obtain 
the ground state energy of the original, or "bare" [8], 
Hamiltonian in the infinite model space with evaluated 
uncertainties. To this end, we incorporate systematic and 
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FIG. 3: (Color online) Calculated ground state energy of ''He 
as a function of Nmax at various values of the oscillator energy, 
hi}, as indicated in the legend. The results are connected 
by straight line segments to guide the eye. The results with 
chiral N3L0 are from Ref. [9]. The results for JISP16 are 
closer to convergence even in rather modest basis spaces. No 
extrapolations are needed in these cases as converged results 
are obtained directly. 



reliable extrapolation tools as needed. 

We compare in Fig. 3 the ground state of "^He using 
the chiral interaction, N3L0 [3, 9], with JISP16 [10, 11] 
and plotted as a function of Nmax to ascertain conver- 
gence rates. All our results include the Coulomb interac- 
tion between protons. Each point represents the ground 
state energy from an A^maa;-truncation of the full infi- 
nite matrix problem. Hence, all points are strict up- 
per bounds on the exact ground state energies for the 
respective Hamiltonian. Fig. 3 shows that JISP16 pro- 
vides a faster convergence rate and a ground state en- 
ergy in closer agreement with the experimental energy of 
—28.296 MeV. JISP16 also produces spectra and other 
observables in light nuclei that are in reasonable accord 
with experiment [10]. Indeed, this interaction was de- 
signed to possess these specific properties while retaining 
an excellent description of the NN data. The interaction 
is non-local but this is no limitation for NCFC which 
also preserves all the symmetries of the underlying NN 
interaction. 

With JISP16 for our NN interaction, we perform ab 
initio calculations of the ground state energies of ^H, ■^H, 
''He, ^He, ^Li, *He, ^^C and i^Q. The three hghtest nuclei 
serve as test cases to demonstrate that the extrapolation 
methods, using results in limited basis spaces, are able 
to predict the fully converged results and to demonstrate 
that our assessed uncertainties are realistic. We limit 
ourselves to examples for which a sufficient set of results 
could be achieved within our current computational re- 
source limits. 



III. FINITE BASIS SPACE EXPANSIONS 

Our results in finite basis spaces satisfy the varia- 
tional principle and show uniform and monotonic con- 
vergence from above to the exact eigenenergy with in- 
creasing Nmax- That is, the results for the energy of the 
lowest state of each spin and parity, at any Nmax trunca- 
tion, are upper bounds on the exact converged answers 
and the convergence is monotonic with increasing Nmax- 
This guarantee of monotonic convergence from above to 
the exact energy facilitates our choice of extrapolating 
function. 

We carefully investigate the dependence of the results 
on the basis space parameters, Nmax and HQ. Our goal 
is to achieve independence of both of these parameters as 
that is a signal for convergence — the result that would 
be obtained from solving the same problem in a complete 
basis. 

Before proceeding, let us explain some additional fea- 
tures of the many-body regulator, Nmax- As introduced 
above, Nmax is the maximum number of oscillator quanta 
shared by all nucleons above the lowest HO configura- 
tion for the chosen nucleus. Its use allows us to factor- 
ize eigenfunctions into intrinsic and center-of-mass (CM) 
components for ease of eliminating spurious CM motion 
effects on all observables. One unit of oscillator quanta is 
one unit of the quantity (2n + I) where n is the principle 
quantum number and I is the angular quantum number. 
If the highest HO single-particle state of this lowest HO 
configuration has Nq HO quanta, then Nmax + Nq iden- 
tifies the highest HO single-particle states that can be 
occupied within this many-body basis. Note that since 
Nmax is the maximum of the total HO quanta above the 
minimal HO configuration, we can have at most one nu- 
cleon in such a highest HO single-particle state. 

The precise method of achieving the factorization of 
the CM and intrinsic components of the many-body 
wavefunction follows a standard approach, sometimes re- 
ferred to as the "Lawson" method [12]. In this method, 
one selects the many-body basis space in the manner de- 
scribed above and adds a Lagrange multiplier term to the 
many-body Hamiltonian X{Hcm — ^fifl) where Hcm is 
the Harmonic Oscillator Hamiltonian for the CM motion. 
With A chosen positive (10 is a typical value), one sep- 
arates the states of lowest CM motion (OSi) from the 
states with excited CM motion by a scale of order Xhfl. 
The resulting low-lying states have wavefunctions that 
are assured to have the desired factorized form. 

It is important to note that our NCFC results for the 
ground state energy for A = 2,3 & 4 are obtained di- 
rectly as we achieve sufficient independence of Nmax and 
HQ. For the other nuclei studied here, we characterize the 
approach to convergence by the dependence of results on 
both Nmax and hfl and investigate the shape of that con- 
vergence in detail. The degree of residual dependence on 
these two parameters provides a measure of the differ- 
ence from the exact result, an estimate of the numerical 
uncertainty in the extrapolation. 
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We employ the parallel-processor code "Many-Fermion 
Dynamics — nuclear" (MFDn) [13] that sets up the 
many-body basis space, evaluates the many-body Hamil- 
tonian matrix, obtains the low-lying eigenvalues and 
eigenvectors using the Lanczos algorithm, and evaluates 
a suite of experimental observables. Working in the 
single-particle HO m-scheme, the lowest 15 states here 
are usually obtained with 300—600 iterations, depending 
on Nmax and the nucleus involved. The required number 
of iterations grows with Nmax- 

The largest matrix we diagonalize for this work corre- 
sponds to ^^O in the Nmax = 8 space with a basis dimen- 
sion about 1 billion. We obtain the lowest 8 eigenstates 
and a suite of observables in 4.5 hours on 12,090 proces- 
sors using the Franklin supercomputer at the National 
Energy Research Supercomputer Center (NERSC). The 
second largest case is ^^C with a basis dimension of about 
600 million for which we obtain the lowest 15 eigenstates 
and a suite of observables in 2.3 hours on 12,720 proces- 
sors using the Jaguar supercomputer at Oak Ridge Na- 
tional Laboratory (ORNL). The above times correspond 
to calculations at a single value of For calculations 
as a function of hfl in the same basis spaces, we use in- 
ternally generated and stored index arrays amounting to 
many terabytes of data so that the second and subse- 
quent hil values each take about 2/3 the time as the first 
case. 

Following the completion of the calculations reported 
here, further speedups have been accomplished with the 
code so that the above-mentioned times are reduced by 
a factor of 2 in future calculations of the same type [14]. 



IV. EXTRAPOLATING THE GROUND STATE 
ENERGY 

A. Simple illustration 

Here we illustrate the convergence properties for the 
nuclear ground state energy in a HO basis with a sim- 
ple model. While the properties of the HO basis are 
useful for many purposes, such as the exact treatment 
of the center of mass motion and the ease of transform- 
ing between relative and single-particle coordinates, the 
asymptotic HO wavefunctions are gaussians while wave- 
functions of finite nuclei will have exponentially decreas- 
ing amplitudes at large distances. Correct long-range 
behavior is important for precision evaluation of energies 
and for many other experimental properties such as elec- 
tromagnetic moments and transition rates. To achieve 
converged long-range observables, we expect to require 
an optimal choice of the Ml value and sufficiently large 
values of Nmax to generate good asymptotic properties. 
To investigate these issues, we evaluate the properties 
of fermions in a finite phenomenological potential, the 
picture that underlies the successful nuclear shell model, 
using a HO basis expansion. 

Consider the properties of a single Slater determinant 



composed of the lowest ^-particle orbits of the Saxon- 
Woods central potential plus a nuclear spin-orbit po- 
tential. This corresponds to a standard mean field de- 
scription of the nucleus with A nucleons and we refer 
to this simple model as the "Extreme Single Particle 
Model" (ESPM). Instead of solving for the single-particle 
states by numerical integration, we diagonalize the 1- 
body model Hamiltonian in a HO basis to simulate the 
procedures of a no-core finite HO basis calculation. 

We adopt ^^C as an example and we perform this di- 
agonalization as a function of Nmax and Hn. In this way, 
we are studying how the lowest s-state and p-state solu- 
tions depend on Nmax and Ml. We identify Nmax with 
the maximum value of the HO quanta (2n -|- I) retained 
in the HO basis expansion so that Nmax = corresponds 
to a pure HO approximation, Nmax = 2 employs 2 basis 
HO functions for the occupied Os-states and Op-states, 
Nmax = 4 employs 3 basis HO functions for the occu- 
pied single-particle states, etc. Note, there is a difference 
between the use of Nmax in this model problem from 
our NCFC approach. Here, since we work entirely in a 
single-particle basis, all particles have simultaneous ac- 
cess to the range of basis states dictated by Nmax while in 
our NCFC, the many-body basis is restricted so that as 
one particle takes more quanta, the remaining particles 
take fewer quanta. 

It should be noted that our model problem more 
closely simulates the traditional "Full Configuration In- 
teraction" (FCI) approach used in quantum chemistry 
where orbits are equally accessible in the many-body ba- 
sis states up to some single-particle cutoff. The origin of 
this difference is our need to retain an exact treatment 
of the center-of-mass motion in our no-core methods for 
finite nuclei (NCFC and NCSM). In spite of this diflFer- 
ence in cutoffs, we find that this simple model is useful for 
illustrating how proper asymptotic wave function prop- 
erties influence convergence rates for self-bound nuclei in 
the HO basis. 

We adopt a Saxon- Woods central, U{r), and spin- 
orbit, Uso{r), potentials similar to a standard choice [15] 
where, for simplicity, we use the same parameters for the 
neutrons and the protons of ^^C 



U{r) 



1 + exp((r 



S-L 



R)/ao) 

2 



Uso 



1 d 

r dr 1 + exp((r — R)/aso) 



with R = r^A^/^ and [y;^^ = 2.0 fm^. Our parame- 
ters selected for this demonstration are [/q = — 32MeV, 
ro = 1.25 fm, oq = 0.65 fm, Uso = 15 MeV and 
o-so ~ 0.47 fm. For the protons we add the Coulomb 
field of a uniformly charged sphere of radius R. 

We add the energies of the occupied orbits, taking into 
account degeneracies, to obtain the total energy of the 
system in the ESPM. The resulting ESPM ground state 
energy for ^^C is displayed in Fig. 4 as a function of h£l for 
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FIG. 4: (Color online) Calcnlated ground state energy of 
in the Extreme Single Particle Model (ESPM) as a function of 
hil and N„iax , the maximum value of the HO states' quanta, 
2n + I, used in expanding the Saxon-Woods single-particle 
states. The curve closest to convergence corresponds to the 
value N,nax = 20 and successively higher curves are obtained 
with Nmax decreased by 2 units for each curve. 




FIG. 5: (Color online) Calculated ground state energy of i^C 
in the Extreme Single Particle Model (ESPM) as a function 
of Nmax for selected values of HQ used in defining the basis 
states. The points correspond to h^l values ranging from 5 
to 30 MeV in 2.5 MeV increments. The curves represent an 
exponential plus a constant fit to each set of results at fixed 
hil, excluding the Nmax = result. Each point carries equal 
weight. 



a range of N^ax values up to and including N„iax = 20. 
The line segments connect the results calculated at se- 
lected values of H^l. Here we observe a pattern that is 
typical of our no-core basis space results presented be- 
low — a sequence of curves with energy decreasing as a 
function of increasing N^ax, consistent with the property 
dictated by the variational principle. With increasing 
Nmax, one approaches convergence signaled by achieving 
simultaneous independence of both Nmax and hfl. In the 
ESPM, we achieve the total energy converged to within 
10 keV at N^ax = 20 over the range /ifJ = 6-11 MeV. 
For the optimal value of h^l = 7 MeV, we achieve a total 
energy to within 170 keV of the exact answer already at 
Nmax = 8, yielding an upper bound within 0.25% of the 
exact result. 

For the purpose of exploring potential extrapolation 
tools, we use the results of the ESPM to map out the 
convergence pattern of the total energy in the present 
work. In subsequent efforts, we will investigate other ob- 
servables in a similar fashion. Thus, we present in Fig. 5 
the total energy as a function of Nmax at fixed values of 
HQ spanning the minima of the curves in Fig. 4. Specifi- 
cally, the points correspond to Ml values ranging from 5 
to 30 MeV in 2.5 MeV increments. We find that, once we 
exclude the Nmax = result, the calculated points ap- 
pear to represent an exponential convergence pattern. To 
confirm this, we fit an exponential plus constant to each 
set of results as a function of Nmax, excluding Nmax ~ 0, 
and the resulting fits are displayed in Fig. 5 as smooth 
curves. That is, for each set of points at fixed hfl, we fit 
the ground state energy with three adjustable parameters 



using the relation 

Egs{Nmax) = aexp{-cNmax) +£'gs(oo). (1) 

In these fits, we assigned equal weight to each point and 
perform a regression analysis. 

Overall, we conclude that the exponential plus con- 
stant fits the results rather well. Thus, one observes that 
the HO basis provides a rapidly converging sequence of 
total energies in the ESPM, one well-represented by expo- 
nential convergence in Nmax towards the the asymptotic 
total energy, Egs{oo). It appears reasonable to expect 
this convergence pattern of the HO basis treatment of 
the ESPM to be representative of HO basis expansion 
behavior in our no-core applications, we will adopt this 
functional form as a foundation for further developing our 
extrapolation methods below. In the following sections, 
we will use additional arguments to improve upon this 
tool and test it in light nuclei where converged results 
are obtained directly. 

We note that a similar exponential behavior for HO 
basis space calculations of a cold trapped Fermi gas has 
been observed [16]. In that case, the same type of single- 
particle-space regulator was employed as we use here in 
the ESPM application. 

The exponential plus constant has also been employed 
as an extrapolation tool in more conventional shell model 
studies as well [17]. In those applications, the variable is 
the matrix dimension rather than Nmax- 
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FIG. 6: (Color online) Calculated ground state energy of 
as a function of the oscillator energy, HQ, for selected values 
of Nmax used in defining the basis states. The curve closest 
to experiment corresponds to the value Nmax ~ 20 and suc- 
cessively higher curves are obtained with Nmax decreased by 
2 units for each curve. The curves are formed by straight line 
segments joining calculated results. 



B. NCFC test case: deuteron 

Next, wc turn to the NCFC calculations for light nuclei 
using JISP16, where we can achieve nearly exact results 
in large model spaces. In this and the following subsec- 
tions, we investigate the convergence rates for the ground 
state energies as a function of N,nax and hfl for ^H, ■^H, 
and ''He. We discuss two extrapolation methods, which 
allow us to obtain estimates of the converged NCFC re- 
sults from finite model spaces. We also introduce the 
assessed uncertainties for our extrapolated results. 

The sequence of curves in Fig. 6 for illustrates the 
trends we encounter in our calculations when evaluating 
the ground state energy with JISP16. Our purpose with 
is only to illustrate convergence trends and test the 
extrapolation tool since the exact answer is also avail- 
able from the direct solution of the Schroedinger equa- 
tion [18], and agrees with experiment. The Nmax — 18 
curve reaches to within 9 keV of this exact result; the 
Nmax = 20 curve reaches to within 5 keV. We note that 
the weak binding of leads to a slow progression of 
the curves towards independence of hfl and contrasts the 
stronger binding situation obtained for '^He shown below 
in Fig. 13. 

1. Global extrapolation method 

We use these results to test our "global extrapola- 
tion method" (referred to as "extrapolation A" ) as illus- 
trated in Fig. 7. Here, we fit only four calculated points 
at each value of hn in the range Nmax = 6—12, represen- 
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FIG. 7: (Color online) Calculated ground state energy of '^H 
as a function of Nmax at hil values that bracket the best upper 
bound as described in the text. The smooth curves are fits by 
Eq. (1) to the 4 data points shown at each value of HQ with 
each point weighted equally. 



tative of the limited results that we expect to encounter 
in slightly heavier systems. We select the values of h^l to 
include in the analysis by first taking the value at which 
the minimum (with respect to hD,) occurs along the high- 
est Nmax curve included in the fit, then taking one hfl 
value lower by 2.5 MeV and three hfl values higher by 
successive increments of 2.5 MeV. Since the minimum oc- 
curs along the Nmax = 12 curve at = 12.5 MeV as 
shown in Fig. 6, this produces the 5 curves spanning a 
range of 10 MeV in hil shown in Fig. 7. We perform a 
linear regression for each sequence of four points at fixed 
hQ using an independent exponential plus constant for 
each sequence, and observe a small spread in the extrap- 
olants as is evident in Fig. 7, which is indicative of the 
uncertainty in this method. 

We recognize that this window of results in fiO values is 
arbitrary. Our only assurance is that it seems to provide 
a consistent set of extrapolations in the nuclei examined 
in the present work. 

For the global extrapolation we chose sets of four 
points due to a desire to minimize the fluctuations arising 
from certain "odd-even" effects already visible in Fig. 4. 
These effects are most pronounced in weakly bound sys- 
tems, and may be attributed to the fact that HO wave- 
functions fall off too fast: wavefunctions of finite nuclei 
decrease exponentially at large distances. In order to 
mimic such an exponential decrease with HO basis func- 
tions, one needs HO basis functions with both even prin- 
ciple quantum number n (even number or nodes in radial 
wavefunction) and with odd principle quantum number 
n (odd number or nodes in radial wavefunction). Since 
N = {2n + 1), a set of four successive Nmax points (with 
Nmax even) implies we incorporate two highest allowed 
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HO single-particle states with even values of the princi- 
pal quantum number and two highest allowed with odd 
values. Thus, a set of four consecutive Nmax points in- 
stead of three points (the minimal number of points for 
an exponential extrapolation) averages out some of these 
"odd-even" effects. We will come back to this point when 
we discuss the extrapolation method B, using sets of 
three points at fixed Ml. 

Next, we consider what weight to assign to each calcu- 
lated {NmaxTMl)-^o\T\t. The fits in Fig. 7 are obtained 
with equal weights for each of the points. However, as 
Nmax increases, we are approaching the exact result from 
above with increasing precision. Hence, the importance 
of results grows with increasing Nmax and this should be 
reflected in the weights assigned to the calculated points 
used in the fitting procedure. With this in mind, we 
adopt the following strategy: define a chisquare function 
to be minimized and assign a (Tjv„,„^ to the gs energy 
at each A^max value that is the change in the calculated 
energy from the next lower TYmax value 



^A^TTiax = {Egs{Nmax) — Egs{Nmax — 2)) 



To complete these sigma assignments, the sigma for the 
lowest A'niax point on the iV,„ax curve is assigned a value 
three times the sigma calculated for the second point on 
the same fixed-fi.51 trajectory. As a final clement to our 
global extrapolation strategy, we invoke the minimization 
principle to argue that all curves of results at fixed hfl 
will approach the same exact answer from above. Thus 
all curves will have a common asymptote and we use that 
condition as a constraint on the chisquare minimization. 

When we use exponential fits constrained to have a 
common asymptote and weights based on the local slope, 
we obtain curves close to those in Fig. 7. The differences 
are difficult to perceive in a graph so we omit presenting 
a separate figure for them in this case. It is noteworthy 
that the equal weighting of the linear regression leads to 
a spread in the extrapolants that is still modest. 

The sequence of asymptotes for the ^H ground state 
energy, obtained with the global extrapolation, by using 
successive sets of four points in Nmax and performing our 
constrained fits to each such set of four points, is shown 
in Fig. 8 as extrapolation A. We employ the independent 
fits such as those in Fig. 7 to define the uncertainty in 
our asymptotes. In particular, we define our uncertainty, 
or estimate of the standard deviation for the constrained 
asymptote, as one-half the total spread in the asymp- 
totes arising from the independent fits with equal weights 
for each point. On rare occasions, we obtain an outlier 
when the linear regression produces a residual less than 
0.999 that we discard from the determination of the to- 
tal spread. Also, on rare occasions, the calculated upper 
uncertainty reaches above the calculated upper bound. 
When this happens, we reduce the upper uncertainty to 
the upper bound as it is a strict limit. 




FIG. 8: (Color online) Extrapolated ground state energies and 
variational upper bounds from each set of four (extrapolation 
A) or three (extrapolation B) successive Nmax values as a 
function of the largest value of Nmax in each set. Error bars 
are dominated by the uncertainties in the extrapolations, and 
are obtained as described in the text. Note the expanded scale 
and the reasonable consistency of the extrapolated results: for 
Nmax > 10 all but one are within their uncertainty range of 
the exact answer. 



One may worry that the resulting extrapolation tool 
contains several arbitrary aspects and we agree with that 
concern. One recourse is to cross-check these choices with 
the solvable NCFC cases in the present subsection and 
following subsections. We seek consistency of the con- 
strained extrapolations as gauged by the uncertainties es- 
timated from the unconstrained extrapolations described 
above. Indeed, our results such as those shown in Fig. 8, 
demonstrate that consistency. The deviation of any spe- 
cific constrained extrapolant from the result at the high- 
est upper limit Nmax appears well characterized by the 
assigned uncertainty. 



2. Extrapolation at fixed HQ 

In addition, we also employ an extrapolation at fixed 
values of hfl using only three successive values of Nmax, 
the minimal number of points for such an extrapolation 
(referred to as extrapolation B). Under the assumption 
that the convergence is indeed exponential, such an ex- 
trapolation should get more accurate as Nmax increases; 
the difference between the extrapolated results from two 
consecutive sets of three Nmax values is used here as our 
estimate of the numerical uncertainty associated with the 
extrapolation. 

In Fig. 9 we illustrate this extrapolation for ^H based 
on calculations with Nmax = 8, 10, and 12. As we can 
see, this extrapolation gives h^l dependent results. We 
therefore consider the value of hfl where the extrapo- 
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FIG. 9: (Color online) Calculated ground state energy of 
for Nmax = 8, 10, 12, and the extrapolated ground state 
energy using method B, as function of the oscillator energy, 
h^. Error bars are obtained from the difference with the 
extrapolation using Nmax = 6, 8, and 10 calculations. 




Oscillator Energy (IVIeV) 

FIG. 10: (Color online) Calculated ground state energy of 
as a function of the oscillator energy, ftfi, for selected values 
of N,nax- The curve closest to experiment corresponds to the 
value Nmax ~ 20 and successively higher curves are obtained 
with Nmax decreased by 2 units for each curve. 



lation is most stable (i.e. for which the difference be- 
tween the extrapolated value and the result at the high- 
est NmCix is minimal) as the best or most reliable Ml for 
this extrapolation method. This hft value is usually at 
or slightly above the variational minimum. 

Since this extrapolation uses sets with only three N„iax 
points, the "odd-even" effects may be significant, in par- 
ticular for weakly bound nuclei. This is indeed what we 
find for as seen in Fig. 8. Nevertheless, within the 
estimated error bars, the results are consistent with ex- 
trapolation method A and with the exact result. In addi- 
tion, as we proceed to applications in heavier nuclei and 
more deeply bound nuclei, this extrapolation becomes 
more stable and useful, as we will see below. 



C. More NCFC test cases: A = 3, 4 

The ground state energies of "^H using JISP16 are 
shown in Fig. 10 as a function of the HO energy for 
the same sequence of basis spaces as for ^H. We again 
observe a converging sequence of upper bounds with an 
indication of a small amount of underbinding compared 
with experiment. Wc note that the curves show a greater 
region of approximate independence of hTl than found in 
the case of as may be expected from the stronger bind- 
ing relative to the first breakup threshold in the present 
situation. 

We use the case of ■^H to illustrate again the utility of 
the global extrapolation (A). Results of independent fits 
with equal weights for each calculated point are shown 
in Fig. 11 to demonstrate the nearly identical asymptote 
when results are available at sufficiently high Nmax val- 
ues. We depict in Fig. 11 both the experimental value 
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FIG. 11: (Color online) Calculated ground state energy of 
as a function of Nmax at values that bracket the best upper 
bound. The smooth curves are fits of an exponential plus a 
constant [see Eq. (1)] to the 4 data points shown at each value 
of h^l. There is a 5 keV spread in these asymptotes that is 
used to assign the uncertainty to the asymptote derived from 
the constrained fit as described in the text. The asymptote of 
the constrained global fit (extrapolation A) is quoted in the 
figure and the experimental result is shown for comparison. 



and the asymptote of the global extrapolation using the 
four point sequence up to Nmax = 18. The span of hfl 
values is selected in the same manner as in the case, 
the procedure we will use throughout this work. 

As for the case of ^H, we also present the sequence 
of extrapolated results for the '^H ground state energy 
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FIG. 12: (Color online) Extrapolated ground state energies 
and variational upper bounds from each set of four (extrap- 
olation A) or three (extrapolation B) successive Nmax values 
as a function of the largest value of Nmax in each set. Error 
bars are dominated by the uncertainties in the extrapolations, 
and obtained as described in the text. Note the consistency 
of the extrapolated resuhs. 



in Fig. 12 using both extrapolation methods A and B, 
together with the variational bound. Both extrapolation 
methods appear to be consistent with each other, and 
give numerical error bars which decrease with increasing 
Nmax- The extrapolation B, using only three successive 
Nmax points at fixed HQ, shows a rather strong "odd- 
even" effect. Nevertheless, all extrapolated results agree, 
within error bars, with each other, and with our best 
results at Nmax = 20. These results are also quoted in 
Table I. We conclude that JISP16 undcrbinds by 
approximately 113 keV. 

Our calculations for '^He show a similar convergence 
pattern as those for ^H. At Nmax = 20, our results are 
within a few keV of full convergence, as can be seen from 
Table I, and we find JISP16 underbinds '^He by about 
52 kcV. 

As a final test of our extrapolations wc consider *Hc. 
We present our calculations as a function of h^l at fixed 
values of Nmax hi Fig. 13. The results clearly indicate 
rapid convergence in both Nmax and hfl; at Nmax = 16 
the ground state energy is converged to within 1 keV 
over the range 20 MeV < hfl < 25 McV. Furthermore, 
the fully converged NCFC ground state energy is within 
3 keV of the experimental energy as shown in Table I. 
As illustrations of our extrapolations, we demonstrate in 
Fig. 14 the independent fits used to assess uncertainties 
of extrapolation A based on Nmax = 2 to Nnax = 8 
results. In addition. Fig. 17 shows extrapolation B at 
fixed values of hil using Nmax = 8, 10, and 12 results. 

We also confirm that the dependence on Nmax at fixed 
Ml is nearly a pure exponential as illustrated best in 
Fig. 15 where we show a wider range of the calculated 




15 20 25 30 
Oscillator Energy (MeV) 



r 

35 



40 



FIG. 13: (Color online) Calculated ground state energy of 
^He as function of the oscillator energy, HQ, for a sequence 
of Nmax values. The curve closest to experiment corresponds 
to the value Nmax = 16 and successively higher curves are 
obtained with Nmax decreased by 2 units for each curve. 
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FIG. 14: (Color online) Calculated ground state energy of 
^He for Nmax = 2—8 at selected values of hQ. Each set of 
4 points is fit, using equal weights, with an exponential plus 
constant [see Eq. (1)] producing the solid curves. Half the 
resulting spread in the asymptotic values is used to deter- 
mine the uncertainty assigned to the first point in Fig. 16 for 
extrapolation A as described in the text. 



results. Here, we provide regression analyses for each 
set of results spanning Nmax = 2—16 at fixed Ml values 
ranging from 15 to 35 MeV. For both basis space pa- 
rameters, this is a significantly wider range of parameter 
values than we use in our applications below. 

We present the NCFC results of both extrapolation 
methods in Fig. 16 along with the experimental and vari- 
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TABLE I: Binding energies in MeV of nine nuclei and of the first excited 0"'' states in C and O from experiment and theory. 
The experimental values are from Refs. [19-24]. The uncertainties in the rightmost digits of an extrapolation are quoted in 
parenthesis. The bounds for the binding energies follow from the variational upper bounds for the ground state energies. The 
rightmost column provides the uppermost value of A^max used in the quoted extrapolations. 



Nucleus (,/^) 


Exp 


Extrap (A) 


Extrap (B) 


variational bound 


Max(iV, 


(1+) 


2.225 


2.223(5) 


2.226(6) 


2.220 


20 


'H(i + ) 


8.482 


8.369(1) 


8.3695(25) 


8.367 


18 


'He (1+) 


7.718 


7.665(1) 


7.668(5) 


7.663 


18 


"He (0+) 


28.296 


28.299(1) 


28.299(1) 


28.298 


18 


"He (0+) 


29.269 


28.68(12) 


28.69(5) 


28.473 


14 


^Li (1+) 


31.995 


31.43(12) 


31.45(5) 


31.185 


14 


*He (0+) 


31.408 


29.74(34) 


30.05(60) 


28.927 


12 


i^C (0+) 


92.162 


93.9(1.1) 


95.1(2.7) 


90.9 


8 


(0+) 


84.508 


80.7(2.3) 






8 


I'^O (0+) 


127.619 


143.5(1.0) 


150 (14) 


134.5 


8 


i«0 (0+) 


121.570 


130.6(7.6) 
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FIG. 15: (Color online) Calculated ground state energy of 
"He for Nmax = 2—16 for JISP16 at selected values of hO,. 
Each set of 8 points at fixed hO, is fit by Eq. (1) producing 
the solid curves. Each point is a true upper bound to the 
exact answer. The asymptotes Egs (oo) are the same to within 
35 keV of their average value and they span the experimental 
ground state energy. 



FIG. 16: (Color online) Extrapolated ground state energies 
and variational upper bounds from each set of four (extrapo- 
lation A) or three (extrapolation B) successive Nmax values as 
a function of the largest value of Nmax in each set. Error bars 
represent the assessed uncertainties in the extrapolations and 
are obtained as described in the text. Note the consistency of 
the extrapolations: the exact answer is well within the uncer- 
tainty range of the extrapolations for all Nmax points, with 
the uncertainty diminishing with increasing Nmax- 



ational upper bound energies. In this case the results pro- 
duce very rapid convergence with uncertainties that drop 
precipitously with increasing Nmax as seen in the figure. 
We note that the error bars conservatively represent the 
extrapolation uncertainties since all the extrapolated re- 
sults are, within their error bars, consistent with each 
other and with the fully converged NCFC result. The 
largest Nmax points define the results quoted in Table I, 
a ground state ovcrbound by 3 ± 1 kcV. 

We have seen in this section that the NCFC results for 



three light nuclei provide sufficiently converged ground 
state energies to allow us to test our extrapolation meth- 
ods and confirm the validity of their assigned uncertain- 
tics. In what follows, we present NCFC calculations for 
five nuclei using both extrapolations A and B. The five 
nuclei selected for this initial application consist of sta- 
ble and unstable even nuclei that span the p-shell. We 
include two loosely-bound nuclei, ^He and *He, antici- 
pating that they will provide challenges for achieving a 
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FIG. 17: (Color online) Calculated ground state energy of 
''He for Nmax = 8, 10, 12, and the extrapolated ground state 
energy using method B, as function of the oscillator energy, 
hfl. Error bars are obtained from the difference with the 
extrapolation using Nmax = 6, 8, and 10 calculations. 



converged ground state energy. 
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FIG. 18: (Color online) Calculated ground state energy of ®He 
as function of the oscillator energy, ftfi, for selected values of 
Nmax- The figure displays also the experimental result, the 
common asymptote from the global extrapolation (A) and 
the experimental threshold for the a + 2n breakup. The curve 
closest to experiment corresponds to the value Nmax ~ 14 and 
successively higher curves are obtained with Nmax decreased 
by 2 units for each curve. 



D. NCFC results for ''He, ^Li, and '^He 

By comparing Figs. 6, 10 and 13 we observe clearly 
the marked correlation between binding energy and con- 
vergence rate — the more deeply bound ground states 
exhibit greater independence of hTl at fixed Nmax- A 
more complete picture of this correlation is seen below, 
for example, with the ^He and ^Li results where one ob- 
serves that the relevant energy scale governing the rate 
of convergence is the binding with respect to the nearest 
threshold. Our physical intuition supports this correla- 
tion since wc know the asymptotic tails of the bound 
state wave functions fall more slowly as one approaches 
a threshold for dissociation. This same intuition tells us 
to expect Coulomb barriers and angular momenta to play 
significant roles in this correlation. 

Consider first the weakly-bound nucleus ^He presented 
in Figs. 18, 19, and 20. In the largest basis spaces 
achieved, Nmax = 14, we obtain net binding with respect 
to the breakup threshold as seen in Fig. 18. However, the 
results appear farther from convergence than the case for 
^He in the same N^ax spaces depicted in Fig. 13. Note 
the proximity of the ^He -I- 2n breakup threshold to the 
calculated ground state energies that suggests the im- 
portance of achieving results at Nmax = 16 as well as 
obtaining improved nuclear Hamiltonians that better re- 
produce the experimental binding of ^Hc. 

To examine this situation in some detail, we present in 
Fig. 19 results as a function of Nmax that covers a range 
of 15 MeV in Ml values. We also present the linear re- 
gression analyses for the range of Nmax = 2—14 to show 




FIG. 19: (Color online) Calculated ground state energy of ''He 
for Nmax ~ 2—14 for JISP16 at selected values of hfl. Each 
set of points at fixed h^l is fit by Eq. (1) using equal weights 
producing the solid curves. Each point is a true upper bound 
to the exact answer. The resulting asymptotes Egs{oo) are 
the same to within 600 keV of their average value. The figure 
displays the experimental result and the common asymptote 
from the global extrapolation (A) as described in the text. 
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FIG. 20: (Color online) Extrapolated ground state energies 
and upper bounds from sets of successive Nmax values as a 
function of the largest value of Nmax in each set. Uncertain- 
ties are determined for each value of fiQ as described in the 
text. Note the consistency of the extrapolated results as they 
fall well within their uncertainty ranges along the paths of 
converging sequences. 



FIG. 21: (Color online) Calculated ground state energy of ®Li 
as function of the oscillator energy, h£l, for selected values of 
Nmax- The curve closest to experiment corresponds to the 
value Nmax ~ 14 and successively higher curves are obtained 
with Nmax decreased by 2 units for each curve. The figure 
displays also the global extrapolation (A) and the threshold 
for the a + d breakup. 



that the exponential fit appears to maintain its validity. 
Thus, we proceed with the extrapolation methods as de- 
veloped and tested in previous sections and display the 
results in Fig. 20. 

Since we use only three Nmax points for our extrapola- 
tion B, at fixed values of Hi}, we also include extrapolated 
results based on Nmax = 2, 4, and 6; with an error es- 
timated based on the difference between the Nmax = 6 
calculation and the extrapolated result. Such an extrapo- 
lation, based on a rather small model space, can be useful 
for larger nuclei, and it turns out to be quite reasonable. 
However, we do see a rather significant "odd-even" effect 
with extrapolation B for this nucleus. 

The results from the highest value of the upper limit in 
Nmax are provided in Table I. We again observe consis- 
tency in the results the the global and fixed Ml extrap- 
olations. The JISP16 interaction yields about 600 keV 
underbinding in ^He which implies the theoretical pro- 
ton RMS radius will likely be too large compared with 
experiment. 

Next, we turn our attention to ^Li and present our 
NCFC results using extrapolations in Figs. 21, 22 and 
23. As in the ^Hc case, there is a low-lying threshold for 
breakup — here about 1.47 MeV above the experimental 
ground state. Both our Nmax = 12 and 14 curves drop 
below this threshold over a range of h^l values as seen in 
Fig. 21. 

Our global extrapolation for ^Li is depicted in Fig. 22 
where we select the case with the highest upper limit in 
Nmax to portray. The minimum in hfl at Nmax = 14 
occurs at hfl = 20 MeV. According to our global extrap- 
olation, we then perform the constrained fit on the results 
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FIG. 22: (Color online) Calculated ground state energy of Li 
for Nmax ~ 8—14 at values of h^l that span the minimum at 
Nmax = 14. Curves define the fits using the global extrapo- 
lation (A) that produces a common constant, the asymptote, 
labelled by "Extrap". 



in the span ;irj = 17.5-27.5 MeV in 2.5 MeV increments 
to obtain the fits shown in Fig. 22. The asymptote, the 
extrapolant (—31.43 MeV), becomes the last data point 
on the right in Fig. 23 where the uncertainty is obtained 
in the manner described above. We also perform our 
extrapolation at fixed values of HQ, and find results con- 
sistent with our global extrapolation, see Fig. 23. 

The extrapolated results are entered in Table I and 
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FIG. 23: (Color online) Upper bounds and extracted asymp- 
totes for the ground state energy of ''Li from each set of suc- 
cessive Nmax values as a function of the largest value of Nmax 
in each set. Uncertainties in the asymptotes are determined 
as described in the text. Note the consistency of the extrapo- 
lated results as they fall well within their uncertainty ranges 
along the paths of converging sequences. 



FIG. 24: (Color online) Calculated ground state energy of ''He 
for Nmax ~ 6, 8, 10, 12, and the extrapolated ground state 
energies using method B, as well as the result from extrap- 
olation method A, with error bars. For comparison, we also 
show the experimental value, and the threshold for a + 4n 
breakup. 



compared with experiment. We find that Li is under- 
bomid by about 560 keV with the JISP16 interaction, 
similar to the amount of underbinging for ^He; the rate 
of convergence and error estimate in our final answer are 
also similar. However, the "odd-even" effect we found 
with extrapolation B for the other nuclei (including ^He 
below) is absent in ^Li. 

Next, let us consider another weakly-bound nucleus, 
®He. In Fig. 24 we show our results for the ^He ground 
state energy as function of HQ for several values of Nmax, 
together with the extrapolated results at fixed hCl and our 
result from the global extrapolation. The extrapolated 
energies and their uncertainties are presented in Fig. 25 
as function of Nmax , together with the variational upper- 
bounds. For comparison, we also show the experimental 
value, and the thresholds for a + An breakup and for 
^He + 2?? breakup. 

Clearly, the results are not as well converged as those 
for lighter nuclei, because we are limited to a smaller 
model space, Nmax = 12. Nevertheless, the variational 
upperbound on the ground state energy is well below the 
a + An threshold. Furthermore, the extrapolations are 
consistent with each other, and the error bars decrease 
with increasing Nmax ■ Our final NCFC result is not only 
below the a +An threshold but also below the experimen- 
tal ^He + 2n breakup threshold, even taking into account 
the uncertainty in the extrapolation, as can be seen from 
Fig. 25 and Table I. Compared to the experiment, we 
find that *^He is underbound by about 1.6 ± 0.4 MeV 
with JISP16. 




FIG. 25: (Color online) Upper bounds and extrapolated 
ground state energies of *He from each set of successive Nmax 
values as a function of the largest value of Nmax in each set . 
Uncertainties in the asymptotes are determined as described 
in the text. 



V. NCFC RESULTS FOR ^^C AND ^'^O 



Having illustrated the application of our methods with 
a range of light nuclei, we finally turn to heavier nuclei, 
and perform NCFC calculations for ^^C and ^^O. For 
these nuclei, we can only go to Nmax = 8, so we have to 
rely on the extrapolation methods. 
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FIG. 26: (Color online) Calculated ground state energy 
as function of the oscillator energy, hil, for selected values of 
Nmax- The curve closest to experiment corresponds to the 
value Nmax = 8 and successively higher curves are obtained 
with Nmax decreased by 2 units for each curve. 



A. Extrapolating ^^C 

The dimension of the model space for ^^C with N„iax = 
8 and limited to total M = states in the m-scheme 
is 594,496,743; and the total number of nonzero matrix 
elements in the lower triangle of the extremely sparse 
many-body Hamiltonian matrix is 539,731,979,351 with 
NN interactions only. Thus, storage of one vector in this 
model space requires 2.4 GB, and storage of the lower 
triangle of the matrix requires 4.3 TB. The dimension 
of the Nmax = 10 basis space is 7,830,355,795, which is 
beyond our present capabilities. 

In Fig. 26 we show our results for ^^C for Nmax = 
through Nmax = 8. Since Nmax = is generally not 
very reliable for our extrapolations, we only have only 
the extrapolation from the Nmax = 2—8 results. In or- 
der to illustrate the details of our uncetainties, we depict 
in Fig. 27 the linear regression analyses of our results 
spanning the minimum in hfl obtained at Nmax = 8. 
Our global extrapolation A produces a ground state en- 
ergy of 93.9 ± 1.1 MeV, whereas the extrapolation B at 
fixed /lil = 27.5 MeV (where it is most stable) gives 
95.1 ± 2.7 MeV. Given the "odd-even" effect which of- 
ten plagues extrapolation B, in combination with the 
smaller error bar obtained with extrapolation A, we con- 
clude that JISP16 produces a binding energy of about 
94 MeV with an uncertainty of 1% to 2%; or in other 
words, it overbinds by about 1.8 MeV. 

For a speculative application, we also consider the first 
excited 0+ state of ^^C, the "Hoylc state" or "triple- 
alpha" state as it has come to be known. Since ex- 
perimentally, this state, with i^Hoyic = —84.51 MeV, is 
just above the threshold for breakup into three alpha's, 
3Ea = —84.89 MeV, it may be poorly converged. On 



FIG. 27: (Color online) Calculated ground state energy ofi^C 
for Nmax = 2—8 at selected values of HQ, as described in the 
text. For each ftSl the data are fit by Eq. (1). These indepen- 
dent asymptotes Egs{oo) provide a measure of our uncertainty 
within the global extrapolation (A). The figure displays the 
experimental ground state energy and the common asymptote 
obtained in the global extrapolation. 
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FIG. 28: (Color online) Calculated energy of ^'^C first excited 
0^ state (Hoyle state) for Nmax = 2—8 at selected values of 
hO, as described in the text. For each HQ the data are fit by Eq. 
(1). These independent asymptotes Egs{oo) provide a mea- 
sure of our uncertainty within the global extrapolation (A). 
The figure displays the experimental energy and the common 
asymptote of the global extrapolation. 



the other hand, our calculations for both ^He and ®He 
at Nmax = 2—8 arc above breakup into alpha plus neu- 
trons, but the extrapolations from these points produce 
results with error bars of about 1 MeV, and agree with 
our best calculations at Nmax = 12. Encouraged by these 
results for ^He and ®He at Nmax = 8, we apply the global 
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FIG. 29: (Color online) Calculated ground state energy of 
as function of the oscillator energy, hfl, for selected values of 
Nmax- The curve closest to experiment corresponds to the 
value Nmax = 8 and successively higher curves are obtained 
with Nmax decreased by 2 units for each curve. 



extrapolation method A to the first excited 0+ state of 
^^C. The calculated results and extrapolation are shown 
in Fig. 28 and summarized in Table I. Our extrapolation 
gives i^Hoyie = —80.7 ± 2.3 MeV, corresponding to an 
excitation energy of 13 ± 3 MeV, compared to an experi- 
mental excitation energy of 7.654 MeV. It remains to be 
seen how reliable the extrapolation is for this (and simi- 
lar) states. One may even expect this extrapolation to be 
unreliable as solutions obtained in our present, very lim- 
ited, basis spaces may not accommodate all the essential 
physics of such excited states. Nevertheless, assuming 
that our error estimates are realistic, our conclusion is 
that JISP16 overbinds the ground state of ^2C by an 
MeV or two, but underbinds the first excited 0+ state by 
about 2 to 6 MeV. When combined, that means it pro- 
duces an excitation energy that is significantly too large. 



B. Extrapolating ^''O 

Finally, we consider -'^^O with Nmax = 8. The dimen- 
sion of the corresponding model space is 996,878,170; and 
the total number of nonzero matrix elements in the lower 
triangle of the many-body matrix is 805,811,591,748 with 
NN interactions only. Thus, storage of one vector in this 
model space requires 4.0 GB, and storage of the lower 
triangle of the matrix requires 6.5 TB. 

The results for ^^O are shown in Figs. 29 and 30, and 
summarized in Table I. As in the case of ^^C, we only 
attain the results through Nmax = 8 with our current 
capabilities. Thus, we have a single extrapolant using 
Nmax = 2, 4, 6, and 8 (extrapolation A), or use extrap- 
olation B for Nmax = 4, 6, and 8. The extrapolation 
B gives a lower extimated ground state energy with a 
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FIG. 30: (Color online) Calculated ground state energy of 
^^0 for Nmax = 2—8 at selected values of HQ, as described 
in the text. For each hi}, the data are fit to an exponential 
plus a constant, the asymptote. The figure displays the ex- 
perimental ground state energy and the common asymptote 
obtained in the global extrapolation (A). 



significantly larger uncertainty (10%) compared to ex- 
trapolation A. We anticipate that this difference is due 
to the "odd-even" effect we have seen in most light nu- 
clei using extrapolation B, and expect that the results 
obtained with extrapolation A are more realistic for this 
case. Also note that the variational upper bound on the 
ground state energy is below the experimental ground 
state. That implies that JISP16 produces an overbind- 
ing of at least 7 MeV, and more likely about 15 to 18 
MeV for ^^O. 

For another speculative application, we also consider 
the first excited 0+ state of ^^O, believed to have a 
significant cluster structure. Experimentally, this state 
is very close to threshold for breakup into ■^'^C plus 
an alpha particle. Applying extrapolation A, we find 
^'excited = —130.6 ± 7.6 McV, showing an even larger ex- 
trapolation uncertainty than for the Hoyle state in ^^C. 
Compared to experiment, we find an excitation energy 
of 13 ± 8 MeV, compared to 6.05 MeV experimentally. 
Given the large uncertainty in the extrapolation for this 
state, we cannot draw any conclusions without results in 
larger model spaces regarding this excited state. 



VI. CONCLUSIONS AND OUTLOOK 

We present in Table I a summary of our results, using 
extrapolations, performed with methods introduced here 
and compare them with the experimental results. In all 
cases, we used the calculated results to the highest Nmax 
available with the bare JISP16 interaction. In the cases of 
the lightest nuclei, the extrapolations were rather modest 
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FIG. 31: (Color online) Calculated first excited 0"'" state en- 
ergy of cis function of the oscillator energy, HQ, for se- 
lected values of Nmax ■ The curve closest to experiment corre- 
sponds to the value N,nax = 8 and successively higher curves 
are obtained with N,nax decreased by 2 units for each curve. 
The figure displays the experimental excited state energy and 
the common asymptote obtained in the global extrapolation 
(A). 



as nearly converged results were obtained directly. The 
uncertainties apply to the least significant digits quoted 
in the table. 

Our overall conclusion is that these results demon- 
strate sufficient convergence is achieved for ground state 
energies of light nuclei to allow extrapolations to the 
infinite basis limit and to estimate their uncertainties. 
Thus, we have achieve ab initio no-core full configuration 



(NCFC) results for these nuclei with our chosen Hamil- 
tonian. 

The convergence rate reflects the short range proper- 
tics of the nuclear Hamiltonian. Fortunately, new renor- 
malization schemes have been developed and applied that 
show promise for providing suitable nuclear Hamiltonians 
with good convergence properties [25]. Additional work 
is needed to develop the corresponding NNN interaction. 
Also, further work is in progress to develop extrapola- 
tion tools for the RMS radii. Of course, the RMS radii 
present a greater challenge since they are more sensi- 
tive than the energies to the asymptotic properties of 
the wavef unctions. 
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